A Double- Transition Scenario for Anomalous Diffusion in Glass-Forming Mixtures 
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We study by molecular dynamics computer simulation a binary soft-sphere mixture that shows 
a pronounced decoupling of the species' long-time dynamics. Anomalous, power-law-like diffusion 
of small particles arises, that can be understood as a precursor of a double-transition scenario, 
combining a glass transition and a separate small-particle localization transition. Switching off small- 
particle excluded-volume constraints slows down, rather than enhances, small-particle transport. 
The data are contrasted with results from the mode-coupling theory of the glass transition. 
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Transport properties in disordered media are impor- 
tant in a wide range of applications from biophysics to 
geosciences. Intriguing behavior arises from 'fast' species 
moving through a dense host system, such as power-law- 
like dynamical conductivities of ion-conducting melts [1] . 
Likewise, 'anomalous diffusion' appears in many amor- 
phous media: mean-squared displacements (MSD) that 
grow like Sr 2 ~ (with some positive (i < 1) over large 
time windows, instead of obeying Einstein's law for or- 
dinary diffusion (/i = 1), are seen in biophysical tracer 
experiments investigating cellular environments [2-4] , in 
zeolites [5, 6], gels [7, 8], amorphous semiconductors and 
photoconductors [9], or specially confined colloidal sus- 
pensions [10-12]. 

These systems can be thought of as mixtures composed 
of a small (fast) species and slow (big) host particles 
providing a highly complex confining structure (called 
'molecular crowding' in biophysical literature). One way 
to deal with this, is to simplify the discussion to stochas- 
tic lattice gases and single tracers moving in a random 
environment [13-15], invoking as a reference point the 
single-file diffusion of non-overtaking particles, Sr 2 ~ i 1 / 2 
[16-18]. Such modeling obviously leaves out two aspects: 
the dynamics leading to a time-scale separation in the 
first place, and interactions among the carrier particles. 

In order to highlight the remarkable features arising 
from dynamical many-body effects in anomalous diffu- 
sion, we investigate a binary, disparate-size soft-sphere 
mixture. We show how anomalous diffusion can be in- 
terpreted as a high-density phenomenon, specifically as 
the approach to a double-glass transition. Many-body 
interactions manifest themselves in a striking way in the 
dynamics of the small species: releasing excluded-volume 
constraints, their mobility is reduced at long times, rather 
than enhanced. 

The appearance of two kinds of glasses - one where 
both particle species freeze, one where the smaller one 
stays mobile - has been predicted [19-21] using mode- 
coupling theory of the glass transition (MCT) [22], and 
indicated in colloidal experiments [23, 24] and molecular- 
dynamics (MD) simulations [25]. MCT qualitatively ex- 



plains fast-ion diffusion in sodium silicate melts [26] as 
a precursor of this scenario. The two transitions have 
different microscopic origins: while the slow dynam- 
ics of the larger species is dominated by caging on the 
nearest-neighbor scale, the single-particle dynamics of 
the smaller species exhibits a continuously diverging lo- 
calization length. This latter leads to the appearance 
of power-law-like anomalous diffusion. A similar transi- 
tion also appears in a recent extension of MCT where big 
particles are immobile from the outset [27, 28]. 

The exemplary off-lattice model for particle localiza- 
tion is the Lorentz gas (LG), a single classical point par- 
ticle moving between randomly distributed, fixed hard- 
sphere obstacles. At a critical obstacle density, the par- 
ticle undergoes a localization transition understood as a 
critical dynamic phenomenon [29]. Close to the transi- 
tion, a power-law asymptote for the mean-squared dis- 
placement is explained by continuum percolation theory, 
as demonstrated in recent extensive simulations [30-32]. 
We shall embark on the subtle connection between the 
LG and true binary mixtures below. 

We performed molecular-dynamics (MD) simulations 
of an equimolar binary mixture of purely repulsive 
soft spheres, with interaction potential V a p(r) = 
4e a/3 [(r / 'a a p)~ 12 - (r / 'cr Q/3 )~ 6 ] +e Q/3 for r < r_ = 1 x ^a af} 
(zero else), a, (3 € (1, s). Diameters are chosen additively, 
ct q /3 = {&aa + app)/2, (Tii setting the unit of length, and 
a ss = 0.35. Nonadditive energetic interactions further 
decouple the species, en = e ss = 1 but e s i = 0.1. We 
set the temperature /cb? 1 = 2/3, and all masses equal, 
mi = m s = 1: no time-scale separation exists between 
the two species at short times. 

The smoothened potential V(r) x [(r — r_)//i] 4 /[l + 
(r — r_)//i] with h = 0.005fJn provides continuity of 
energy and forces at the cutoff r_ . Newton's equations 
of motion were integrated for N\ — N s = 1000 parti- 
cles with the velocity form of the Verlet algorithm (time 
step St = 0.005/y4s in units t = [miafjen] 1 / 2 ). To 
avoid crystallization, big-particle diameters were sampled 
cquidistantly from the interval an £ [0.85,1.15], retain- 
ing <7si = (1 + <7 ss )/2. At each number density p, four 



2 




10 

1.2 1.6 2.0 2.4 2.8 3.2 

P°n 

FIG. 1: MD-simulated self-diffusion coefficients for small (As) 
and large particles (Di) in a disparate-sized binary mixture. 
The solid line is a power-law fit tx (p c — p) 1 , where p c — 2.23 
and 7 = 2.1. Inset: partial static structure factors S a /3(q) 
at two different densities as indicated, as functions of qp~ 1//3 . 
The dotted line is S ss (q) for non- interacting small particles at 
p = 2.19. 



independent runs were performed. Up to p < 2.296(7^ , 
the system was fully equilibrated, requiring equilibration 
runs over at least 10 6 and up to 2 x 10 8 time steps, fol- 
lowed by production runs of the same length. During 
equilibration, temperature was held constant by coupling 
the system periodically to a stochastic heat bath; produc- 
tion runs were done in the microcanonical ensemble. At 
the highest density p = 4.215 erjj 3 , over 10 9 time steps 
were performed. No runs showed signs of demixing or 
equilibrium phase transitions. 

Figure 1 displays the self-diffusion constants D a ob- 
tained from the simulated mean-squared displacement 
(MSD), Sri 0) = - ^(0)) 2 > for a singled-out par- 

ticle at r^(t) via the Einstein relation, 8r^(t — + oo) ~ 
6D a t, where possible. The diffusion coefficients show a 
decoupling between the motion of large and small par- 
ticles which becomes more pronounced with increasing 
density, due to a faster slowing down in D\ than in D s . 
At p — 2.296, D s is about 2.5 orders of magnitude higher 
than D\, and at p > 2.568, big-particle diffusion has 
ceased over the entire simulation time window. Yet, the 
small-particle MSD still retains a diffusive regime, allow- 
ing us to extract D s > up to p = 3.257. Also shown 
in Fig. 1 is a fit to D\ by a power law as predicted by 
MCT, D ~ D Q (p c - p)i. The fit yields p c = 2.23 for the 
critical density and 7 = 2.1 for the critical exponent. A 
similar fit to D s is not satisfying and yields much larger 7 
and p c - D s (p) also clearly differs from an Arrhenius law, 
although a fit with two exponentials would be possible. 

The slowing down visible in Fig. 1 and discussed in 
the following is purely dynamic; no essential changes 
in the static structure of the system were observed, de- 
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FIG. 2: Mean-squared displacements (MSD), 5r 2 a (t) for the 
large (dotted) and small (solid lines) particles in the simulated 
binary mixture. Densities shown are pofj = 2.0, 2.296, and 
4.215 for a = 1 and pal = 2.0, 2.296, 2.654, 2.837, 3.257, 
3.627, 3.906, and 4.215 for a = s. Symbols mark the crossover 
to ordinary diffusion, where 5r^(t) ~ t (slope indicated by the 
dash-dotted line). 



spite the drastic compression employed. This is demon- 
strated by the partial static structure factors S a f3(q) = 
(J2jk ex PM? ' {r a ,j - rp,k)]), showing little change with 
density if plotted as functions of q* = qp^ 1 ^ 3 to eliminate 
a trivial change in length scale (inset of Fig. 1). 

The MSD corresponding to these transport coefficients 
are shown in Fig. 2. For the big particles, we observe 
a standard glass-transition scenario: a two-step process 
gives rise to a plateau over an increasingly large time 
window, crossing over to diffusion at increasingly large 
time, and at a length scale associated with dynamic 
nearest-neighbor cageing, typically at about 10% of a 
particle diameter (Lindemann's criterion). Indeed, from 
the plateau of the p = 2.296 curve one reads off the 
corresponding cage localization length rf = y/ Srf /6 ~ 
0.06 en, which decreases at larger p due to compres- 
sion. Although this could technically be called anomalous 
diffusion, we reserve that term for the behavior shown 
by the small particles: Around p c , the small-particle 
MSD behave quite differently, with no sign of a two- 
step glassy dynamics. Instead, they show subdiffusive 
growth and cross over to ordinary diffusion at increas- 
ingly large length and time scales when increasing p. 
This indicates that nearest-neighbour cageing is not the 
dominant mechanism for their slowing down. The sub- 
diffusive regime can be described by power-law variation, 
Srf(t) cx i M with some effective < p < 1 that seems to 
decrease with increasing density. 

The small-particle dynamics qualitatively agrees with 
previous MD results [25]. It also agrees with the dynam- 
ics found in the Lorentz gas [30-32]. There, subdiffu- 
sive growth with apparent density-dependent exponents 
p is due to the approach to an asymptotic power law, 
5r^(t) ~ t x , that extends to t — > 00 at the localization 
critical point. Careful simulations [30] could establish 
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FIG. 3: (a) Small-particle MSD with (solid lines) and 
without (dashed) interactions between small particles, den- 
sities as indicated. (b) Logarithmic derivative fi(t) — 
d[\o g 8r 2 s (t)]/d(\ogt). 

x = 2/6.25 for the LG. To estimate a critical exponent 
x from Fig. 2 is tempting, but preasymptotic corrections 
render it impossible. It appears that a description of our 
data in terms of the LG asymptote is not convincing. 

Our binary mixture differs from the LG inter alia 
through the finite density of interacting small particles. 
To establish the effect of this distinction, we switch off 
interactions among small particles, setting e ss = while 
keeping their number constant. Within simulation ac- 
curacy, structure and dynamics of the big particles are 
unchanged in this 'transparent-small' mixture. 

Figure 3 compares the small-particle MSD of the two 
systems. Initially, the transparent small particles show 
weaker localization, intuitively expected as they have 
larger free volume available. This trend prevails at low 
densities. Surprisingly, at high p, switching off inter- 
actions leads to significantly slower diffusion compared 
with the fully interacting case as t — > oo. 

This is emphasized by the lower panel of Fig. 3: the ef- 
fective exponent fi(t) = d[\ogSr^(t)]/d(\ogt) crosses over 
from //(0) = 2 (ballistic short-time motion) to //(oo) = 1 
for ordinary diffusion or //(oo) = for arrested parti- 
cles. For the LG model, /x(t) x for increasingly large 
time windows close to the localization treshold. No clear 
plateaus are seen in our data, but switching off small- 
small interactions at fixed density p systematically re- 
duces fi(t) at long times. For comparison, we have in- 
dicated in Fig. 3 the predictions x = 1/2 for single-file 
diffusion and x = 2/6.25 = 0.32 from the LG model. 

One could rationalize this finding as follows: excluded- 



FIG. 4: (a) Mean-squared displacements 5r^(t) obtained 
from mode-coupling theory with MD-simulated S(q). Solid 
(dashed) lines: small-particle MSD with (without) interac- 
tions between small particles, densities as labeled. Dotted: 
big-particle MSD for the highest density shown, (b) Corre- 
sponding logarithmic derivative fJ,(t), as in Fig. 3. 



volume constraints dominate the exploration of the 
small-particles' local surroundings. In the long run, how- 
ever, the exploration of all the cul-de-sacs for the nonin- 
teracting small particles in the frozen structure becomes 
vastly less effective than an interaction-mediated trans- 
port: interacting small particles have a larger probabil- 
ity to visit spaces where another small particle has just 
left, thereby channeling motion along 'preferential paths' 
known from ion-conducting melts [33]. 

No exact results are known for our binary mixture. 
MCT describes the slow glassy dynamics well, but a 
small-wave-vector divergence forbids its application to 
the localization transition [34, 35], and hence to a dis- 
cussion of x. Still, the predicted double-transition sce- 
nario [19-21] is in line with our results. Thus, ignoring 
the asymptotic behavior close to the localization tran- 
sition, let us focus on the effect of making the small 
particles 'transparent' in the binary mixture. The MCT 
equations for Sr^(t) are completely determined once the 
S a f)(q,t) are known; we take them from MD (cf. inset of 
Fig. 1). The difference in interactions enters the theory 
only through these quantities. To focus on effects arising 
from finite wave numbers, we only show MCT solutions 
on a finite grid qi — (i ~ l/2)Aq<7f 1 with i = 1, . . . 120, 
Aq = 0.2 and additional \ow-q cutoff qo = 4<Tj _1 . Dis- 
cretization of q in fact yields imct = 1/2 [21, 36]. 

Figure 4 shows the Sr^(t) for the two binary mix- 
tures considered. Remarkably, the theory reproduces 
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three qualitative trends seen in the MD data, Figs. 2 
and 3: (i) while the big-particle MSD exhibits the ordi- 
nary glassy two-step behavior (localization length around 
O.lci), Sr^(t) shows a different signature in the time win- 
dows accessible to the simulation. This is the precursor 
of the double-transition scenario, (ii) At low densities, 
transparent-small diffusion (dashed) is faster than the 
one for interacting small particles (solid lines), as ex- 
pected from the reduced scattering frequency. This also 
holds for high densities at intermediate times, (iii) For 
large t and high p, transport of transparent small par- 
ticles is much slower than in the fully interacting case. 
Within the MCT picture, the latter arises from a small 
shift of p c w 2.293 (7 w 2.9) to lower p, rendering trans- 
port slower at fixed p. From our MD data, we cannot 
rule out whether this MCT picture is qualitatively cor- 
rect, or whether one may indeed surmise that the change 
in dynamics is due to a cross-over from single-particle dy- 
namics akin to the Lorentz-gas model to a many-particle 
interaction- assisted transport. 

Let us summarize the main results. We studied a 
disparate-size mixture of purely repulsive soft spheres 
whose dynamics can be understood as the approach to 
two distinct, purely dynamical arrest transitions: (i) 
an ordinary glass transition connected with big-particle 
transport, where small-particle diffusion does not van- 
ish, and (ii) a localization transition for small-particle 
transport at a higher density. As a precursor, a window 
over increasingly large length scales appears in the small- 
particle mean-squared displacement, exhibiting power- 
law anomalous diffusion, 6r^(t) cx t x . 

This naturally explains an order-of-magnitude decou- 
pling between diffusion coefficients, rendering our binary 
soft-sphere mixture a minimal model for fast ion trans- 
port in amorphous materials. Further experiments on 
specifically catered colloidal suspensions [23, 24] would 
be highly promising. 

The anomalous diffusion in our binary mixture is a 
many-particle phenomenon: upon switching off interac- 
tions between the small particles, effective power-law ex- 
ponents appear to decrease. As a consequence, excluded- 
volume interactions between the small particles acceler- 
ate their transport in the binary mixture. This is remark- 
able, since in the high-density regime one usually expects 
excluded volume to hinder individual particle motion; to 
understand the origin of this effect remains as a challenge 
for future work. 

We acknowledge discussions with T. Franosch and 
F. Honing and thoughtful comments by W. Gotze. ThV 
is funded through the Helmholtz-Gemeinschaft, young 
investigators group VH-NG 406. 
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